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Diagrammatic approach proposed many years ago for strong correlated Hubbard model is devel- 
oped for analyzing of the thermodynamic potential properties. The new exact relation between such 
renormalized quantities as thermodynamic potential, one-particle propagator and correlation func- 
tion is established. This relation contains additional integration of the one-particle propagator by the 
auxiliary constant. The vacuum skeleton diagrams constructed from irreducible Green's functions 
and tunneling propagator lines are determined and special functional is introduced. The properties 
of such functional are investigated and its relation to the thermodynamic potential is established. 
The stationary properties of this functional with respect to first order changing of the correlation 
function is demonstrated and as a consequence the stationary properties of the thermodynamic 
potential is proved. 

PACS numbers: 71.27.+a, 71.10.Fd 



I. INTRODUCTION 

The Hubbard model is one of the most important mod- 
els for the electron of solids which describes quantum 
mechanical hopping of electron between lattice sites and 
their short ranged repulsive Coulomb interaction. 

This model was discussed by Hubbard M j n order to 
describe a narrow-band system of transition metals and 
has been revised to investigate the properties of highly 
correlated electron systems such as the copper oxide su- 
perconductors and others. 

Hubbard model exhibits various phenomena includ- 
ing metal-insulator transition, antiferromagnetism, ferro- 
magnetism and superconductivity. This model assumes 
that each atom of the crystal lattice has only one electron 
orbit and the corresponding orbital state is non degener- 
ate. 

The Hamiltonian of Hubbard model is a sum of the 
two terms 

H = H° + H', (1) 

where H° is the atomic contribution, which contains the 
Coulomb interaction term U and local electron energy e 
on the atom 

i 

Hh = V en^> + Un^r , (2) 
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€ = 1 — /i, = Ci> C—t> 
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and hopping Hamiltonian 

* 3 

t(t-t) = t*(f-t), t(0) = 0. 



Here Ci> ) are the creation (annihilation) electron 

i a 1 " _^ 

operators with local site i and spin a. Because in the 

thermodynamic perturbation theory we shall use thermal 

averages in a grand canonical ensemble we have added to 

the Hamiltonian (1) the term —/J,N e 

^ = E"f W 
' - ' L a 

i a 

where /i is the chemical potential and N e electron number 
operator. The quantities U and N e are the fundamental 
parameters of the model and because of large value of 
the Coulomb repulsion it is taken into account in zero 
approximation of our theory. The operator H' , which 
describes hopping of the electrons between sites of the 
crystal lattice is regarded as a perturbation. 

For investigating this model new physical and mathe- 
matical concepts and techniques have been elaborated. 
We can only enumerate some of them. There are 
many analytical approximations as Hubbard approxi- 
mation, noncrossing approximation (NCA), slave-boson 
method, Dynamical Mean Field Theory (DMFT), com- 
posite and others projection operator methods and ap- 
proaches. Also numerical simulations of thermodynamic 
quantities and density of states have been performed 
by various methods. Some exact results for one- and 
two-dimensional Hubbard model are known. Each ap- 
proximate method has its advantages and disadvan- 
tages. A short and comprehensive reviews of the meth- 
ods can be found in papers and books I 2 " 7 !. Besides 
these approaches we must enumerate also the special di- 
agram techniques elaborated for strongly correlated elec- 
tron systems. Stasyuk I s ! , Zaitsev I 9 1, Izyumov I 10 ' and 
coworkers have developed a diagram technique for Hub- 
bard model, based on the disentangling of nondiagonal 
Hubbard operator x mn ou t of a time-ordered products 
of other such operators. Due to a more complicate al- 
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gebra of Hubbard transfer operator than that for Fermi 
operators the essential features of this technique remain 
till now poorly developed. 

The other diagrammatic approach around the atomic 
limit also has been proposed for Hubbard model both in 
normal I 11 * 12 ! and in superconducting state 1 13 1. This the- 
ory introduces the Generalized Wick Theorem (GWT) 
that uses cumulant expansion of the statistical average 
values for the products of the Fermion operators. The 
GWT takes into account the fact that Hamiltonian H° 
is nonquadratic in fermion operators due to Coulomb in- 
teraction. This last circumstance is responsible for the 
appearance of the nonvanishing site cumulants called ir- 
reducible Green's functions. These new Green's functions 
take into account all the spin, charge and pairing fluctua- 
tions of the system. The perturbation formalism around 
atomic limit has the advantage that local (atomic) phys- 
ical properties can be evaluated exactly and that the 
transfer of this local information to neighboring sites due 
to the kinetic mobility of the conduction electrons can be 
handled perturbatively in powers of hopping integral. 

A GWT for chronological averages of products of elec- 
tron operators was formulated later also by Metzner I 14 l . 
Mctzner did not derive a Dyson-type equation for the 
renormalized one-particle Green's function and the role 
of our correlation function has not established. 

In spite of the similarity to our works, the diagrams 
in his approach are quite different. The n - particle cu- 
mulant is represented a 2n - valent-point vertex with at- 
tached n entering and n leaving lines, where as our such 
a cumulant or irreducible Green function is represented 
by a rectangle with 2n vertices with n entering and n 
leaving arrows. The rectangle contours the vertices with 
the same site index but different time and spin labels. 
In addition Metzner investigated the limit of high lattice 
dimensions. 

In this paper we shall develop the diagrammatic the- 
ory proposed before for Hubbard model I 11,12 ! with the 
aim to demonstrate the existence of the relation between 
renormalized quantities of thermodynamic potential and 
one-particle Green's function and also to prove the sta- 
tionary properties of this potential. 

Such theorem was proved firstly by Luttinger and 
Ward I 15 l for uncorrelated systems by using the diagram- 
matic technique of weak coupling field theory. 

The strong coupling diagram theory used by us needs 
new conceptions and new equations and they are used 
to prove stationary property of thermodynamic potential 
for strongly correlated systems. Such a proof has been 
already achieved for Anderson impurity model in paper 

[16] _ 

The paper is organized in the following way. In sec- 
tion II we develop the diagrammatic theory in the strong 
coupling limit and the skeleton diagrams are introduced. 
In Section III we prove the stationary theorem for renor- 
malized thermodynamic potential. Section IV has the 
conclusions. 



II. PERTURBATION TREATMENT 

We shall use the definition of the one-particle Matsub- 
ara Green's functions in interaction representation as in 
paper I 11 * 12 !; 

G(x\x') = (TC^>)CV v (r')t/(/3))o, (5) 

where x stands for ( x , a, r) and index c for (...)q means 
the connected part of the diagrams which appear in the 
right-hand of part of (5). We use the series expansion 
for the evolution operator U{f3) with some generalization 
because we introduce the auxiliary constant A and use 
XH' instead of H': 

P 

U X {0) = Texp(-A j H'{r)dT). (6) 
o 

In the presence of this constant we shall use index A for 
all dynamical quantities as G\(x\x') and so on. At the 
last stage of the calculations this constant will be put 
equal to one. 

In zero order approximation the one-particle Green's 
function G^(a;|a;') is local 

G^{x\x') = S^S aa ,G^(T-T r ) (7) 

and the Fourier representation of the last function is 

P 

G<p{iu n ) = J G^(r)e^dr 
o 

1 e -PEo _|_ e -flE a e -pE w + e -f3E 2 

Zq ^ iui n + Eo - E a + iuj n + Ea — E 2 ' 

(8) 

Z Q = e-t> E ° + e-^ + + e^ 2 , 

E = 0, E a = E w = e, E 2 = 2e + U, 
uj n = (2n+l)n/p. 

Here Eo,E a and E 2 are the energies of atomic sites. 

As has been proved in papers I 11 ' 12 ] propagator (5) has 
the diagrammatic representation depicted on the Fig. 1. 
The irreducible Green's functions of order n are depicted 
with rectangles with 2n vertices. The arrows which enter 
in the vertex point depict annihilation electrons and that 
which go out the created electrons. 

In paper I 11 ! we have introduced the notion of correla- 
tion function Z\(x\x') which is the sum of strongly con- 
nected diagrams containing irreducible Green's functions 
and related to a more convenient function A\(x\x') = 
G°(x\x') + Z x (x\x'). 

In Fig. 1 the fourth and seventh diagrams of the right- 
hand site belong to correlation function. Due to the fact 
that irreducible functions are local and tunneling matrix 
elements have the property t{lfc — ~x~) = in Fig. 1 
are omitted all the diagrams which contain self-locked 
tunneling elements. 
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FIG. 1: First three orders of perturbation theory for one-particle propagator. The thin solid line depicts zero order propagator, 
thin dashed line depicts the tunneling matrix element. The rectangles depict irreducible two-particles Green's functions G% ' lr 
and points are the vertices of diagrams. 



As is seen from Fig. 1 the process of propagator renor- 
malization is accompanied by the analogous process for 
tunneling matrix elements renormalization and replacing 
of the instant quantity Xt(x — x') = Xt(x — X )S(t — t' — 
+ )) by dynamical one T\(x\x') equal to 



T x (x\x') 



\t(-£ - -£')8(t - t' - 0+) 
^\t(x-l)G x {l\2)\t{2-x') 



(9) 
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which in Fourier representation 

t{~£) = —^2e(t)exp(-iltlt), 



K 

■-rf 



x exp [—i k ("af — it 1 ) — ioj n (T — r')] 



has a form: 



T\(t\iuJn) = \Tx(lc\iu) n ) 



Xe{t)(l + Xe(t)G x (t\iLJ n )) (10) 



The renormalized tunneling matrix element T\ really is 
tunneling Green's function and will be depicted as double 
dashed line. T\ is represented by such double dashed line 
multiplied by A. 

Now we introduce the skeleton diagrams which con- 
tain only irreducible Green's functions and simple dashed 
lines without any renormalization. In such skeleton di- 
agrams the thin dashed lines are replaced by double 
dashed lines with realizing the complete renormalization 
of dynamical quantities. 

The skeleton diagrams for correlation Aa function are 
depicted on Fig. 2. 

As can bee seen from Fig. 2 the skeleton diagrams for 
A(fc) function are of two kinds. The first four diagrams 



of this Figure are local and therefore their Fourier rep- 
resentation don't depend on momentum. The last dia- 
gram and other ones with more number of rectangles are 
non localized and their Fourier representation depends 
on momentum. Only the first category of diagram are 
taken into account in Dynamical Mean Field Theory. 

As was proved in papers I 11 ' 12 ! the knowledge of this 
function permits to formulate the following Dyson-type 
equation for one-particle Green's function 



G x (k) 



A A (fc) 



1 - Ae(¥)A A (fc) 



(11) 



,-rr 



Here k stands for ( k , iu) n ) with odd Matsubara frequen- 
cies. Equation (10) and (11) gives us the results 



Tx(k) - 
Tx{k) = 



XTx(k), 



1 - Xe(k)Ax(k) 



(12) 



Equation (12) has a form of Dyson equation for tunnelind 
Green's function, and the role of mass operator Ea is car- 
ried out by correlation function multiplied by auxiliary 
constant A: 



E A (fc) - XAx(k). 



(13) 



In Hubbard I approximation we neglect the correlation 
function Zx(k) and consider function Ax(k) equal to zero 
order Green's function G ■ (k). In this approximation we 
have 



G{(k) = 



G°(fc) 



l-Ae(k)G°(k)' 



which describes two Hubbard energy subbands distanced 
between them for U ^ and without the possibility to 
describe the Mott-Hubbard transition. 
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FIG. 2: The skeleton diagrams for correlation function Aa. Double dashed lines depict the tunneling Green's functions T\ and 
the rectangles depict the irreducible Green's functions. 



III. THERMODYNAMIC POTENTIAL 
DIAGRAMS 

The thermodynamic potential of the system is deter- 
mined by the connected part of the mean value of the 
evolution operator I 11 ' 12 ! 



F 



Let us consider a more general quantity first 



F(X) = Fo-^(Ux(J3))l. 



(14) 



(15) 



and put then A = 1. 

By using the perturbation theory we have obtained the 
first orders of diagrams for (U\((3))q, depicted in Fig. 3. 

In order to obtain the better understanding of these 
diagrammatic contributions we examine the expression 



G x (x\x')Xt(t' - -&)5(t -t' - 0+)6„ 



(16) 



where double repeated indices suppose summation and 
integration. Consequently (16) is equal to 



= -A e( k )G\ a { k \iuj n ) exp(^„0 + ). 



(17) 



k u 



Here we have carried out the integration by time. 

From diagrammatic point of view equation (16) implies 
the procedure of locking of the external lines of propaga- 
tors G\ diagrams depicted on the Fig. 1 with the tun- 
neling matrix element t{Ht' — x) and obtaining in such 
a way the diagrams without external lines similar with 
ones for (U\(/3))q depicted on Fig. 3. These two series of 
diagrams differ by coefficients in front of them. 



In expression (17) the coefficients — before each dia- 
gram are absent, where n is the order of perturbation 
theory. These coefficients are present in Fig. 3. In order 
to restore these ~ coefficients in (17) and obtain the co- 
incidence with (J7 A (/3))o series it is enough to integrate 
by A the expression (17) and obtain: 



0+).(18) 



The expression (18) displayed in a diagrammatic repre- 
sentation coincides exactly with mean value of the evo- 
lution operator: 



E 



d\'Gya{^ ~^'\ -0 + ). (19) 







In Fourier representation we have 

A 



(Ux(m 



dX' 



E 

k ffu. 



c(k)Gx>a(k\iui n ) exp(iw n H 



(20) 



From (15) and (20) we obtain 

A 

F(X) = F - 



1 



k a 



,7> 



x Gx>a{ k \iijJ n ) exp(«w„0 + ). 



(21) 



Using the definition (14), the equation (21) can be writ- 
ten in the form: 



F(X) 



A 



k a 



x T,x'(k) exp(iw n + ). 



(22) 
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T 2 



From (22) we have 



F(A) 



FIG. 3: The first four orders of perturbation theory for {U\(/3))q. 

By using these definitions we can prove the equations 



dX p 
k a 



E a E r A (fc)S A (fc) exp^„0+) 



Tr(T\E\). 



(23) 



In order to have a full system of equations we add to (22) 
the definition of the chemical potential of the system 

N e = E ££G a (^|tw n )exp(«j n + ), (24) 

"I' w« 

where iV e is the electron number. 

Equation (21) establishes the relation between thermo- 
dynamic potential and renormalized one-particle propa- 
gator. This last quantity depends on the auxiliary pa- 
rameter A and (21) contains an additional integration 
over it and is awkward because that. 

We shall obtain more convenient equation for thermo- 
dynamic potential without such integration by A. For 
that we shall introduce a special functional 

Y(X) = Y 1 (\) + Y'(\), (25) 

where 

Y ^ = ~\ E M<#)AA A (fc)-l) 

K ,<j,uj n 

+ T x (k)XA x (k)]exp(iuj n O + ) 1 (26) 

and Y'(X) is constructed from skeleton diagrams without 
external lines and is depicted on Fig. 4 

The dependence on A in the functional Y'(X) is twofold, 
just through the dependence of the renormalized Green's 
functions G x , T A , A A and £ A and through the explicit fac- 
tors A ra in front of each diagram of Y'(X). 



ST x (k) 
6/3Y'(X) 



= —XA\(k) = — £ A (fc), 
= XA x (k) =£ A (fc), (27) 



5T x (k) 

and as a result we obtain the stationary property 
Sf3Y(X) 



ST x (k) 



= 0. 



(28) 



By using the definition (13) of the mass operator T, x (k) 
we can rewrite the functional Yi(A) in the form 

Yi(X) = -~££[ln(e(t)£ A (fc)-l) 



+ T A (fc)E A (fc)]exp( lW „0+), (29) 
and prove the second form of stationary property 
6Y(X) 



5Z x (k) 



= 0. 



(30) 



if we take into account the Dyson equation for tunneling 
function T x (k). 

Now we shall discuss the derivative by A of the func- 
tional Y(X). As was mentioned above there is twofold 
dependence of A through S A (fc) and explicit through the 
factors A™ before skeleton diagrams of Y'(X) Fig. 4. 

Due to the stationary property (30) we obtain 



dY(X) 
dX 



= E 



6Y(X) SX x (k) dY(X) 
5Y, x (k) SX dX 1 



(31) 



dY'(X) , 
dX 1 



Here we have taken into account the equation (30) and 
that that Yi(A) has not explicit dependence of A. 



FIG. 4: The simplest skeleton diagrams for functional Y' (A). Double dashed lines are tunneling functions T\(k). 



By using Fig. 4 for Y'(X) and the definition of A A of 
Fig. 2 it is easy to establish the property: 



r)Y'(\\ 1 
A ^A k = ^££T A (fc)AA A (fc)exp^, i O+) 

k a 

= ^Tr(T\A\). (32) 
Therefore from (31) and (32) we obtain: 



= irr(T A E A ). 

From equation (23) and (33) we have 

dF(X) _ dY(X) 
X ^X~~ X ~dX~' 

and as the consequence we obtain 

F(X) = Y(X) + const. 



(33) 



(34) 



(35) 



Because for A = perturbation is absent and -F'(O) 
F o ,Y(0) = we have 



F(X) = Y(X)+F . 
Now we can put A = 1 and to obtain 

F(l) =Y(1) + F . 
with the stationary property 
SF 



as 



o. 



(36) 
(37) 

(38) 



IV. CONCLUSIONS 



We have developed the diagrammatic theory proposed 
for Hubbard model many years ago and introduced the 



notion of renormalized tunneling Green's function T in 
addition to known before. We have defined correlation 
function A, the mass operator E for tunneling function 
and establish the Dyson equation for them. The mass op- 
erator was found (13) to be equal to correlation function 
for A = 1. 

We have obtained the diagrammatic representation of 
the correlation function through the skeleton diagrams 
which contain the many-particle irreducible Green's func- 
tions G 2 n with all possible values of n— order of per- 
turbation theory and the renormalized tunneling Green's 
functions. 

We have established the relation between renormalized 
quantities of thermodynamic potential and one-particle 
Green's function. This last function has additional de- 
pendence on auxiliary constant A and has to be inte- 
grated over it. 

We have proved the possibility to avoid such integra- 
tion by A and to introduce the special functional Y(X) 
constructed from skeleton diagrams. 

At the final part of the paper we have proved the sta- 
tionary property of this functional and found its relation 
to stationary of thermodynamic potential with respect 
variation of mass operator or full tunneling function. 
This theorem is the generalization of known Luttingcr 
and Word theorem I 15 l proved for weakly correlated sys- 
tems on the case of strongly correlated systems described 
with the Hubbard model. 
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